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Abstract 

In the recent paper (referred to as I), the resonating mean- field theory (Res-MFT) 
has been applied and shown to effectively describe the so-called two-gap supercon- 
ductivity (SC). In I, a state with large quantum fluctuations has been approximated 
by superposition of two MF wave functions composed of Hartree-Bogoliubov (HB) 
wave functions with different correlation structures. Particularly in I, using a suitable 
chemical potential, at T = the two-gap SC in MgB2 has been well described by the 
Res-HBT. Furthermore the Res-HB ground state generated with HB wave functions 
has almost explained value of the ground-state correlation energy in all the correlation 
regimes including an intermediate coupling regime. In the present paper we will apply 
the Res-HBT to the two-gap SC with high critical temperature T c . We will aim at 
constructing a theoretical foundation for phenomenological theories of the two-gap SC 
at T = and finite temperature. In the single-gap case we will find a new formula 
leading to a higher T c than the usual HB's. 
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1 Introduction 



A topical two-gap superconductivity (SC) with critical temperature T c = 39K has been 
recently discovered in MgB 2 [2]. It may be expected to open a new area in the vigorous 
pursuit by the radical spirit of the resonating mean-field theories (Res-MFTs) [3, 4] to develop 
a theoretical framework appropriate to explore the problem of high T c in superconductors. 
In particular, fermion systems with large quantum fluctuations show serious difficulties in 
many-body problems at finite temperature. To approach such problems, Fukutome has 
developed the resonating Hartree-Fock theory (Res-HFT) [3] and Fukutome and one of the 
present authors (S.N.) have extended it directly to the resonating Hartree-Bogoliubov theory 
(Res-HBT) to include pair correlations [4, 5]. In the recent paper (referred to as I) [6], the 
Res-HBT has been applied and demonstrated to effectively describe one of an exciting topics 
of current interest, the two-gap SC. An appearance of Thermal Gap Equation in the Res-HBT 
is a manifestation of the analogy of the Res-HBT with the usual BCS and HBT [7, 8]. 

Before the discovery of high-T c superconductor, much effort had been devoted to raising 
T c of the usual BCS superconductor in the weak coupling regime [7, 8, 9] and to obtaining 
the Eliashberg's critical temperature in the strong coupling regime [10, 11, 12]. The T c for 
MgB 2 is 39K, which is close to or even higher than the upper theoretical value predicted by 
the BCS theory [13]. Even if SC in MgB 2 is phonon- mediated, a model beyond the simple 
BCS model or the Eliashberg model is required. The existence of two energy-gaps in MgB 2 
at T = has been predicted phenomenologically by Kortus et al. [14] and Liu et al. [15] 
employing a BCS-like weak-coupling theory, using the effective a and n two-band model. 
They have obtained A CT = 7.4 [meV] and A 7r = 2.4 [meV]. On the other hand, employing the 
Eliashberg's strong-coupling theory [10], Choi et al. [16] have obtained A CT = 6.8 [meV] and 
A 7r = 1.8 [meV]. Beyond such theoretical great successes, the Res-HB ground state generated 
with HB wave functions which are equivalent to the coherent state representations (CS reps) 
[17], is expected to almost explain the value of the ground-state correlation energy in all the 
correlation regimes including an intermediate coupling regime. 

To go beyond such phenomenologies, we develop a tentative temperature-dependent 
Res-HBT. A temperature dependent variation is made to satisfy diagonalization conditions 
for thermal Res-Fock-Bogoliubov (Res-FB) operators along a way different from the usual 
thermal-BCS theory [18, 19, 20]. We derive formulas for determining T c and the behaviour 
of the gap near T = and T c . In the particular case of an equal magnitude of gaps but 
with two different phases, we find new formulas boosting up T c to a higher value than the 
usual HB's value. Finally we stress to attempt a projection-method approach to a rigurous 
thermal Res-HBT. A partition function in an SO(2N) (TV: Number of fermion states) group 
can be calculated using the projection method proposed in the Res-MFTs [3, 4] 

In Section 2, we give a diagonalization condition for thermal Res-FB operators which is 
the genralization of the condition to the Res-HBT's proposed by the Ozaki's method [21]. 
The diagonalized thermal Res-HB density matrix is expressed in the form of a Fermi-Dirac 
distribution with Res-HB eigenvalues. These lead to a self-consistent thermal Res-HB gap 
equation, from which we get formulas to determine T c and the behaviour of gaps near T = 
and T c . Finally in Section 3, we give a summary and further perspectives. In Appendices, we 
give a resonating mean-field free energy, variation of the resonating mean-field free energy, 
diagonalization of thermal HB density matrix and necessary integral formulas to calculate 
new formulas for the gap at intermediate temperature. 
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2 Thermal resonating HB equation 



The Res-HB eigenvalue equation J 7 , I can be extended to the thermal Res-HB 

eigenvalue equation. We obtain the thermal Res-HB equation in Appendix A. The Res-HBT 
implies that every HB eigenfunction in an HB resonating state has its own orbital-energy. 
We can derive the thermal Res-HB eigenvalue equation which is given in Appendix B. The 
thermal HB density matrix is determined as W rrp [F rp \ = [l + exp{f3J r rp }}~ 1 which is proved 
also in Appendix B. Using a Bogoliubov transformation g\ p and g 2 p in I and (B.9) and (B.10), 
H /r iip[J r i p ] and W^ppSp] in I are diagonalized as follows: 



Wrp 9rp^^rrp\3~ rp\9rp~ 



W 



rp 







l-w. 



rp 



1—w 



rp " 



1 + e 



-fie 



rp 



= 



1,2), 
1 



k R T, 



(2.1) 



(2.2) 



which are the generalizations of the Ozaki's results [21] to the Res-MFT. By making the 
Bogoliubov transformation g rp , eigenvalues e rp are obtained by diagonalization of the ther- 



mal Res-FB operators T rp with additional terms (H[W rrp \ — E)\c, 



rp | 



r : 



1,2) (B.10). The 

thermal HB interstate density matrix in the whole Res-HB subspace is given as the direct 
sum: 



W p [JF P ] = g p W p g\ = @ 2 r=l W rrp [Trp] , W rr [T r ] = g rp W rp g 



t 

rp ' 



(2.3) 



Suppose a tilde thermal Res-HB density operator W\ p for equal-gaps to be 



W lp = 



wl p -h 







wi P -h 



w 



W). 

lp ■ 



W lp J-2 











(2.4) 

Here I2 is the two-dimensional unit matrix. Performing the unitary transformation by g[ p , 

— t(-l) 

we obtain the following thermal Res-HB density matrix W lp . 





ta). 
i P ■ 



^ta)fj>ta)^ta)t 

yip vv ip yip 



9ip 







■w lp )-l 2 



~TU)t 
9i P 



§{l-cos# lp (1-2^)} I 2 {-(+)} x |sin^ lp e-* (l-2wl^)l 2 
_{-(+)} x isin 6 lp e^ (l-2wl^}l 2 ±{l-Kx>s 9 lp (l-2w]W)}l 2 



(2.5) 



The W 2p has the same form as (2.5). The Res-FB operator P^) P f° r s P m -up state with 
upper and lower signs, corresponding to Case I (4.5) and Case II (4.6) in I, is expressed as 

He p -h {+(-)} xFlph 

.{+(-)} xFlp-h -^e p -h _ 



P - 



(2.6) 



P 



,-{ £p +2(H[W]-E^) 



sin -f (cos -| 



COS 6 r 



A 2 

=F [deU 12 ] 



1 



l±[deUi 2 ]5 



(2.7) 
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^=^-iA{^(0)y.arcsi„hQ ± [det,,^}. sp i- jI , 
for Case I (upper sign) and Case II (lower sign) where 



x- 



D 



[det Z12] 2 =exp 



-2X( ())ri^ L) \ ln(l+x 2 )+2a;-arctan(^ 



(2.8) 



(2.9) 



At finite temperature, using the formulas (2.1) and (2.2) we require orrespondence relations 
cos 9 p t =>• cos 0i p and sin 9 pT =>■ sin 9\ v given through 



cos r , 



■pX _i_T-t 



sin 9. 



A T 



24 

e lp 



(l-2iDQ : 



(1-2^) , 



Notice the multiplication factor l—2w\^\ The ej p (= e|) is the quasi-particle (QP) energy: 



At 
> — > — =- 



Fljl-2w\ p ) 



(2.10) 



p 



J r _l £pT + J~l £pT j Equations in L.H.S. of (2.10) are unified into one equation 

in R.H.S. It is easily shown that equation (2.10) plays the role of self-consistency condition 
at T = 0. Dividing numerator and denominator, respectively by (e p + A 2 -,) 3 / 2 , equation in 
R.H.S. of (2.10) is rewritten as 



{el+A T ) 



2J% 

At 



(l-2tft) 



, w'—w 



(4 £p t+^V)(i-24)' 



i P - 



(2.11) 



(^ + A|)f 

Now we demand a new condition for thermal gap equation 



Ep 

which leads to 



e 2 p +A r )2 



2^ 2 



{el+A\ 



■(l-2«$)=0, (2.12) 



e;Rcs(±) 
flUJ 



-N(0)V- arcsinh =p [det 2i 2 ]|}EpA> 
+^ g R r c T s(±) ^ D E p 5 pT A 2 • [det z 12 ]| E P C P = 0, 

— =±2M0) fiWDxl ■ arcsinh f — ){ 2-N(0)V- arcsinh f — )) [ det ^ 12 ^ i 



(2.13) 



Here we calculate the term (H [W] -E™ s ) in (2.7) using; the solutions for the Res-HB CI 
equation obtained in I and denote the result as E^,^. We also give the following definitions 
for E P A>, Ep5 p and EpC p : 

_ 2 



Ep^pj Ep-^ P 5 EpCp — Ep 

Rearranging (2.13), it is cast to 



1 



1 



[el+A T )r el+A T > (ej+A T )l 



(1- 



(2.14) 
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N(0)V 



: arcsinh 



arcsinh (—) At-AtT B„ 1 
l±2iV(0) - } XtJ T T2 - p 7 -[detz 12 ]| 



x 



l±[detzi2 
arcsinh (£) Ar-A T E p Bp 



-,-1 



[det z 12 



(2.15) 



which reduces to the Res-HB gap equation (4.10) in I as T— >0. Using a variable e = £A T 
instead of e, the summations J2 p A p , J2 p B p and EpCp near T = can be computed to be 



Ep^p f X T 



e 



2N(0) Jo (e + 1)^ 



; arcsinh 



iV(0) 



" — 2arctan( — ), 

\x T . 



1 \ Ep Cp 



1 



2iV(0) 



(2.16) 



Substituting (2.16) into (2.15) and near T = approximating as 



\x T . 



7Y 



arctan — ~ — — xt, 



arcsinh 



l+xl 



:1-x xt, (0<rr <l) 



■) ^ arcsinh (— ^ - — (x T -x ), [det Zi 2 ]t — [det z 12 ]l =Q , 
\x / Xq 



(2.17) 



then, near T = we have the gaps for Case I (4.5) in I as 



A^ = A h- 



1 - [det z 12 



2 

T=0 



N(0)V 



,1 + arcsinh 



X 



,A 

i 

[det z 12 }t =0 



1 - l ~ l mow h ° + ' N ^ A "- 

N i°) V l + [det; 12 ]f,„ 



1 + arcsinh 



'hou 



D 



A 



(2.18) 



and Case II (4.6) in I with the aid of [det Zi 2 ]t= — l—27rN (0)hou dx (0 < x 1) easily derived 
from Taylor expansion of (2.9), as 



A£: 



D ' 



Tr-l-f -N(0)V-^-N(0)V) ■ [det z 12 ] T l 



— + A arcsinh I 



iu D \ 



(2.19) 



In the opposite limit T— >T^ (T c for Case I) the gap becomes very small, [det zi 2 ]^ 2 — >1, 
then Tl T ->-AiV(0)Uarcsinh(^a; D /A)/4 and ( J r J £pT +J rt £pT )/2 -> e p /4 if we use the last 
relation in (5.9) and the second one in (5.16) in I. We have an approximate QP energy 

~ ^Je 2 p + {-AN{0)Varcsmh{hu D /A)} 2 /4. This is because the two HB WFs have different 
correlation structures ip 2 = iT and ipi=0. In such a case, returning to the original form of 
the BCS gap equation but with the modified QP energy e^ 1 , the thermal gap equation is 

expressed as l = V/2J2 p (l — 2tD^ /e^ 1 and leads to the integral form 



1= ^ E 4 tanh (8& 



-AN(0)V 
5 



de- tanh 

e 



Sk B T l J ' 



(2.20) 



Introduce a dimensionless variable ylp=e /SksT 1 and its upper- value ^^^d/S/cbT*. Inte- 
grating R.H.S. of (2.20) by parts, it is approximated as follows: 



1 



47V(0)V 



~ In y l Tl -J dy In y sech 2 y = In y]^ + In [——J 



=in (£tl) " in (D ' :Debye temperature ) ' 

where we have used the formula in the textbook [22]. Number 7 is the Euler's constant 
(7~0.5772) and e 7 ~ 1.781. Finally a small rearrangement yields 

^ = 0.283^"^, (2.22) 

which should be compared with the Eliashberg's formula [23] and the usual HB's one for T c 

T c = 1.13O0 D e"W. (2.23) 

The new formula (2.22) gives a high critical temperature, e.g., T C I = 72.87K for N(0)V = 0.25 
and 9d = 700K. This T\ is in contrast to T c obtained by the usual HB's (2.23), i.e., T c = 14.49K 
for the same values of N(0)V and 9 D . 

We are now in a stage to discuss the behaviour of the gap near T c . In the above the modi- 
fied QP energy e—\J e 2 + {— A/V(0) Varcsinh(/iWD/ A)} 2 / 4 plays a crucial role to boost the T c in 
(2.22) comparing with the numerical result in (2.23). Notice the exsistence of the numerical 
factor 1/4. Now let us consider A^ near T\. Using this form of QP energy, the gap equation 
is roughly rewritten as 



cfe-tanh — — — ) ~ / cfe-tanh 



AN(0)V A Jo e XZk B Tj Jo e m B T, 



(2.24) 



from which we obtain 

1 cr 



AN(0)V 



\k R Tir Tl 8vr 2SW VeV ^ T} ^ 



I' ' ( ' (■ \ ' c 

I 1 1 -J. I ) \ ' / I \ ^ - 

X 



k B Tl) ^Vzt-^Q-)}', 

where huD/ksT = e 7 /27r • UiOD/ksT. For details see Appendix C. Using arcsinh(l/a; T ) ~ 
m(2/x T )~-(:r T -2)/2 + •• • (2/x T >l/2), (2.21) and (2.25), we get A^ near T\ as 



A^/X^U Tl ~ T " 



7C(3)iV(0)UV n )\ T\ ■ 



T}-T . . 

(2.26) 



The JT^ — T dependence of A^ is more complicated than the usual one [19, 20]. 
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For Case II, (J-^ £ T + J^ e T ), jF^ t and become infinite simultaneously in the limit 

1 /9 

A T — >0 (xr— >0) due to the existence of 1 — [detz^r in denominator. Then mathemati- 
cal handling for such a problem is too difficult and therefore we can not easily get a for- 
mula for T c n in an analytical way as we did in Case I. Denote T c for Case II as T c n . At 
T~T c n , A^J almost vanishes and 1 - [det z 12 ]t 2 -^2irN(0)hw D x T . Using 7"2 T in ( 2 - 8 )> we 
reach to the following asymptotic forms: — > — (47rAf(0))~W(0)^arcsinh(/iw D /Ar) and 

{P +£pT + P £pT )/2 -> (4^(0))-%/ A T . The QP energy ej (= 

is approximately calculated to be e^ n ~ (47rA^(0)) _1 £ p /A T (0<A T ^C 1). Here we discard 

the contribution from JFX T comparing with the one from (J^ e T + J r l £pT )/2. As was done 
previously, returning again to the original form of the BCS gap equation but with another 
modified QP energy e£ n , the thermal gap equation is obtained as 1 = V/2J2 P (l — 
which also leads to the integral form 



r 

■ N(0)V 
Jo 



e p \2k B T ll -A7rN(0)A^J 

(2.27) 

■ 4^V(0)Ag / g \ 

d£ £ tanh V2^rn.47r7V(0)A«J ■ 

Introduce a dimensionless variable = e/ (2kBT-4irN(0)A^fj and its upper-value = 

huD/(2kBT u -4TrN(0)Alfj. Integrating the last equation in (2.27) by parts, it is approxi- 
mately calculated as 



1 1 , II ! 



which reads 



/ 1 e~< hu D \ ( 1 9 D \ 
n {irN(0)AU2irk B T 11 ) ^7^(0) A" f u ) ' 



^ f 1 1 0r> /! 1 



(2.28) 



A r = ^77^7f ex Pf ^77^777^ n) « ^77^77 I 1 ~ ^77^777^ nj > ( 2 -29) 



7rA"(0)T n 1 4/V(0)VW(0)A^j %N(0)T lI \ 4N(0)VttN (0)A 
from which we obtain an equation to determine A^ very near T c as 



ii2 1 9d a h , 1 #d I 



AM 2 _^A"+^— : = 0. (2.30) 

T 7rN(0)T n T 7rN(0)T n 4N(0)VW(0) V ; 



Then we have a solution for A" as 

A 11 — 1 ° D 1 (on) 

T 7iN{0)f u 4N{0)VirN{0Y 1 ; 

in which at A^ = T c n , the A^ vanishes. Then finally we can determine the critical tempera- 
ture A" c for Case II as 

Tf = — d D N($)V=\M4d D N(Q)V. (2.32) 

7T 



The simple formula (2.32) gives a high critical temperature, e.g., T c n = 198K for N(0)V = 
0.25 and #d = 700K. Finally near T c n can be approximately obtained as 



AH, 



i e D T— t 11 



(2.33) 



2tt vriV(O) T c n T c n ' 

which is linearly dependent on T— T c n . It is very interesting that we could find such a 

dependence of A 1 ^ , comparing with the usual dependence \J T— T c n of Aj.. 

In intermediate temperature region the modified QP energy e is approximated as 

-l 



e<±) = y/^ + A^ |2 ^l±[det zi 2 ]fJJ . When e> A T , each term £ p A p , E P # P and £ p C p in 
(2.14) is approximately computed as 



2N(0) 

AtE.5; 



£ 2 +A 



iV(0) 



2iV(0) 



/ ^ ^tanh — — =— C 3 fy r a; T J , 

! ( £ 2 +A 2 ) 5 V^T; tt 2 V > 



(2.34) 



whose details are presented in Appendix C. Taking only a leading term, A p and C p terms in 
(2.34) are approximated to be 



i 



e 



2N(0) 

A t^2p C p _ 7 C(3) ^,(±) ~j 2 _ 



l±[detzi 2 ]| T y 



2iV(0) 
(±)_21C(3) /e 7 



<u°; / (±)~ y 
--^r(y ( T V) 



2 / X \ 2 

A T 



27f2 V 71 l±[detz 12 ]|/ v^Ty 



, A T = A T iV(0)\/arcsinh| 



(2.35) 



Substituting these results into (2.15), we have 



\ n l±[det z 12 



A x^(0)yarcsinhQj|v(0)Varcsinh(— ^-^l-[det z l2 ] 



2tt 2 



-In I 



7 D 



21C(3) \n i ±[detZl2 ]| T)\e D ) 
2tt 2 / 



-Wo)Varcsinh(— 



21C(3) 



2 /e 7 



(1 [del ::•> , ) 

2 







D 



7Y 



l±[det ^i 2 ]| 



T 



2 T / v?d 



(2.36) 



3 W 



l±[det z 12 ] J 



— x^[det z 12 \t = 0, 
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to be solved analytically for a given T, which is rewritten as 

2 



H x|(iV(0)Varcsinhf— ")} 

* l±[det z 12 ]*/ 1 W J 

iV(0)Varcsinh^-(l-[det zi 2 ]f)} 



21C(3) \n i ± [det^ 12 ]| 
1 



2 /e 7 



- x|[det 2i 2 ]j. « 0, (0<x t <1) 



" 3 l7r l±[detz 12 ] 

from which we obtain an equation to determine A T for a given T as 



7V(0)y- 



1 



— x^arcsinh — 
71 l±[detz 12 ]| Kx t 



1 



2tt 2 



In 



21C(3)^j \n 1± [ detZl2 ]i 



-In 



7 D 



T 



This is classified into the following two cases: 
Case I: [det z 12 ]| « 0.3 



0.436iV(0)\/a; T arcsinh (-^) = - v / 0J82^n(0.436) - In (J^j ■ (^J , 

from which, using the approximate relation arcsinh(l/xT) ~ (xt~ 2)/2, finally 
solution for x T (iV(0)F = 0.25 and ^ = 700K) as 



XT—i- 



V0.436x0.25 



0.436 x 0.25-2vfT782y/ln(0.436) -In 



T 



T 



v700/ V7ooy' 

Case II: 1 — [detzi 2 ]| « 2nN(0)hu D x T and In ^— j arcsinh^— ^ 

(0.283) 2 {A^(0)\/} 2 arcsinh 2 (-^) -0.782 {vriV(0)^ D } 2 ^ arcsinh(-^) 
-0.782{nN(0)hu D } 
whose solution is easily otained as 

arcsinhQ = (0.283)~ 2 {iV(0)\/r 2 {7riV(0)^ D }^ 



\n( y 0.283)-\n{2nN(0)huj D }-\n 



T 



'D 



T 



'D 



0.391{7riV(0)^ D }(^ 



-{(0.391) 2 {nN(0)hco D y^—j 



rp «. Z 

' ] + (0.283 ) 2 {.V(())\ 7S2 



x 



|ln(0.283)-ln{27riV(0)^ Z) }-ln 



i -i 

2 



from which, using again the approximate relation arcsinh(l/:c r ) ~ (xt~ 2)/2, finally we have 
a solution for x T (iV(O)F = 0.25, N(0)hu D = 0M and D = 7OOK) as 



0.391 {3.14x0.01} 



x T = 2 + 2(0.283)- 2 (0.25)- 2 {3. 14x0.01}^ 
- j(0.391) 2 {3.i4x0.01} 2 (^^ + (0.283) 2 (0.25) 2 x0.782 

x |ln(0.283)-ln{2x3.14x0.01}-ln 



T 

700. 



(2.43) 



700. 



It is very interesting to investigate behaviour of the temperature dependence of the gap 
A T for Case I and Case II. 
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3 Summary and further perspectives 



We have concentrated on derivation of thermal gap equations within the framework of 
Res-HBA. From the Res-FB operators T\ and Ti with equal-gaps, we have found the diag- 
onalization conditions for them, which are essentially of the same form as the former one. 
It leads to the self-consistent thermal Res-HB gap equation and makes possible to derive 
the new formulas to determine T c and the gaps near T = and T c . The formula for Case 
I gives a high T], e.g., T\ = 72.87K for N(0)V = 0.25 and Debye temperature 8 D = 700K. 
This is in contrast with T c of the usual HB formula giving T c = 14.49K for the same values 
of N(0)V and 6d- The formula for Case II gives also a high T c n , e.g., T c n = 198K for the 
same values of N(0)V and 9 D . The temperature dependence of the gap near T = and T c 
becomes more complicated than that of the usual HB and Abrikosov descriptions [19, 20]. At 
intermediate temperature, we have got the solutions of At for Case I and Case II. We have 

1 /2 

taken [det zu] T ~0.3 (Case I) and N(0)huD = 0-01 (Case II) to get real solutions. We have 
got, however, a; T = 0.65 (Case I) and 1.05 (Case II) for T = 70K which are a little bit large 
compared with the real solutions. Improvement of such results should be made. Further, it 
would have been better to draw numerical aspects of the temperature dependence in all the 
correlation regimes. This is possible in the near future. 

For unequal two-gaps, it is also possible to realize the above-mentioned diagonalization 
condition for Res-FB operators T rp (r = l,2). Transforming by a unitary matrix g rp , T rp is 
easily diagonalized. Noticing the same correspondence as the one in (2.10), cos ^ rj3 =^cos 6 rp 
and sin^ r p=^sin^ rp , we assume each diagonalization condition (2.11) holds even in this case. 
Then we obtain coupled equations through a function of A 1T and A 2 t expressed as 



4 ( 




(4+ A 't) 3/ H 


A rTj 



w 



t = L (3-1) 



which reduces to equation in R.H.S. of (2.10) if A 1T = A 2T . The quantities J-^ A and J-l± ep 
are given by the equations similar to (5.9) in I but with more complicated forms of A^r and 
A 2 t- For the time being, as was done in the previous section we here also use the function 
(e p + A 2 rT ) 3 / 2 by which we divide numerator and denominator, respectively, in (3.1). After 
equating the numerator to the denominator and using the relation 1— 2wJ p = tanh(e r /2/csT), 
we sum up over p, namely integrate both sides of the equation over e, to achieve the optimized 
conditions. We obtain coupled thermal Res-HB gap equations and reach our ultimate goal 
of computing temperature-dependent two-gaps. Along such a strategy and method, we will 
make a numerical analysis to demonstrate the behaviour of temperature-dependent two-gaps. 

To solve such a problem, we must provide a rigorous thermal Res-HBA. We have an 
expression for partition function in an SOQN) CS rep \g) [17], Tr(e~ l3H ) = 2 N ~ 1 f (g\e~ l3H \g)dg 
(Jalg is the group integration on group SO(2N)). Following Fukutome [3], introducing the 
projection operator P to the Res-HB subspace, the partition function in the Res-HB subspace 
is computed as Tr(Pe _/3i/ ). This kind of trace formula is calculated within the Res-HB 
subspace by using the Laplace transform of e~^ H and the projection method which leads 
to an infinite matrix continued fraction. A thermal variation of the Res-HB free energy is 
carried out after the inverse Laplace transform of matrix elements of e~^ H . This is made 
parallel to the usual thermal BCS theory, which will be given in a separate paper. 
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Appendix 



A Resonating mean-field free energy F[Z] 



Res 



Suppose a matrix Z to be a usual FB operator, introduce a quadratic HB Hamiltonian 
and the usual HB free energy [21] 



H[Z] 



\[c\ c] Z 



c, 



F[Z] = Tr(H> + \^{W [Z] \n(W [Z])}, W [Z] = 



-m\z\ 



T±( e -f}H[z\y ^ 



(A.l) 



which leads to 



F[Z] = (H- H[Z]) Z - - lnTr(e-^^), 
^g) _ TV(e-^//[Z]) 



Tr(e-^l) ' 



Tr(e-^^l) 



(A.2) 



We have another form for this free energy, i.e., a well-known formula expressed in terms of 
a usual HB density matrix W[Z\ as 

= (H) z +~Tr {W[Z} \nW[Z] + (1 2N - W[Z\) ln(l 2iV - W[Z})} , 

T 1 \ ( A - 3 ) 

W[Z] = 

l-K[Z]* l N -R[Z\\ 
and the trace formulas for the pair operators 

Tr \w[Z] (V Q + ±fy a ) | = P Q/3 [Z], Tr [Z\E Pa } = X Q/3 [Z]. (A.4) 

We give 2 by a direct sum oi J- r , Z — J2r=i ®3~ r and assume each J>(r = 1, • • • , n) is a Res- 
FB operator. Then, instead of the above H[Z] in (A.l), we introduce a quadratic Res-HB 
Hamiltonian 



H[Z} Res = - [c\ c, • • • , c\ c, • • • , c\ c] 







T 

J n 





c 




c], 




c 




c], 




c 







(A.5) 



Along the same way as the one in (A.l), the Res-HB free energy can also be defined. 

Now let us introduce a projection operator P (P 2 = P = P<) to the Res-HB subspace, 
P\ty) = \^/ Rcs ) as P = J2rs=i \9r){S~ l ) rs {9s\, where the |g r )'s are HB wave functions and 
S — (S rs ) (= [det Zrs] 1 / 2 ^ is an n x n matrix composed of the overlap integrals and = S. 
Using the projection operator P, we give the Res-HB free energy in the form 



F[Z] Rcs = Tr(W [Z] Rcs H) + -Tr \w[Z} Rcs \n(W [Z] Rcs ) \ , W [Z} Res = 



p e -f3H[Z] Res p 

Tr[Pe~^ W ^Rc~ 



(A.6) 



in which, by making Taylor expansion of lnP = ln{l — (1 — P)} and using P 2 = P, we have 
-^Tr {w[Z} R ^ln(w[Z} Res )} = -(H[Z} R ^ z . Res - - ]nTr(Pe^ H ^). (A.7) 
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A natural extension of the HB free energy to the Res-HB free energy is easily made. Then, 
the Res-HB free energy is given as 



F\Z\ 



Res 



(H - H[Z] Rcs ) Z ; Rcs - i lnTr(Pe-^Res), 



_Tr(Pe-^RosH[Z} Rcs ) _ Tr (P e -^^H) 

{H [A\ Res ) Z - Res = r>„-BH\Z]v \ ' \ H /-Z;Res = 



Tr (p e -/3ff[Z] Re 



Tr(Pe-P H WRe S ] 



(A.8) 
(A.9) 



, n) in which 



Consider the whole Res-HB subspace |# Res ( fc )) = Y% =1 <& k) \g t ), (k = 1, 
the Res-state with index k — 1 and the Res-states with indices k — 2,- • -,n stand for the 
Res-ground one and the Res-excited ones, respectively. In order to determine the thermal 
|<7 r )'s and thermal mixing coefficients c^'s by the variational method, we use a temperature 
dependent Lagrangian with the Lagrange multiplier term E to secure the normalization 
condition (^< k )\^< k )) = 1; 

L[Z\H = ELi E?^=i [^ rs [Z] Rcs ] - • [det ^ 4 fc )*4 fc ). (A.10) 

The variation of (A. 10) is made in a quite parallel manner to the one in the previous works 
[3, 4]. We omit detaied derivations for such parallel cases. From the variation of LfiT]™ 
with respect to c^* for any k, we get a thermal Res-HB CI equation to determine 



£?=i {H [W rs [Z] Rcs ] - £«} ■ [det z rs }^ = 0. 



(ATI) 



We define D rs and D rs as D rs = u s z~l5u\ and D rs = Su s z~}ul, respectively, where z rs = u\u s 
and uj = [fej ,aj] (a r and 6 r :r-th HB amplitudes). The variations of thermal HB interstate 
density matrix H^i^Res and thermal overlap integral [detz] 1 / 2 are given by 

SW rs [Z] Res = D rs ( 1 2N - W rs [Z] Res ) + ( 1 27V - W TS [Z] Res )D rs , ' 



5 [det z r 



D r 



[det z r 



(A.12) 



We obtain also the variation of Hamiltonian matrix element iJ[Wj. s [iT]R es ] as 
^[H/ rs [Z] Rcs ]=iTr{^[H/ rs [Z]R CS ]5H/ rs [Z] Rcs }, 



J[W rs [Z} Res 



Frs;a/3 [Z\ 



Res 



Frs [Z] Res D rs [Z\ R es 
-D* r [Z] Res - F* r [Z] Res _ 

5H[W rs [Z] Res 



5R 



rs;/3a[ 



Res 



D 



rs;a/3 



Z] Rcs = 



5H[W rs [Z] 



Res J 



h a /3 + [a/3\'y5]R rs; s~ f [Z] Rcs , 



= --[ui\f35}K rs . &1 [Z] Res . 



(A.13) 



SK* sr .^[Z] Res 2 

Following I and the Res-HF theory [3], writing L[Z]%£ = ££ =1 E" s =i £™ {k) [£]R CS 4 fc) *cf) 
and £™( k )[^B J e B = {^[W„[2l R e B ]-S (fc) }-[det ^ rs ] 1/2 , from the variation of Lg£ for any k, we 
obtain a thermal Res-HB equation to determine the thermal mean field wave function u r as 



ELiE^i/C«[z] Res 4-cv 



/C«[Z] Res ^(W-^Z] Res )^^^ 

From (A. 14) we can derive a thermal Res-HB eigenvalue equation. See next Appendix. 



>(A.14) 
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B Variation of resonating mean-field free energy F[Z 

The thermal Res-HB coupled eigenvalue equations is expressed as follows: 
[J- r [Z] ResMr ]* = e ri u rh e ri = e rs - £Li {H[W rr [Z} Rcs ] - £«} |4 fc )| 2 , 

T r [z] Rcs =T r [z]l cs = T[w rr [z} Rcs ] ELi I4 fc) l 2 



HB 

Res 



(B.l) 



We call the 2Nx2N matrix J>[i?] Res the thermal Res-FB operator. From now let us de- 
note Wrs^Res, J>[Z] Res and K,^[Z] Rcs simply as W rs , T T and /C^, respectively. First, due 
to idempotent-like product properties W rs W rr =W rs and W sr W rr =W rr , we have important 
relations Kf}W rr =K^ and K^W„={H[W sr }-E^) W„- [det z sr ]i Next multiplication of 
the second equation in (B.l) by W rr from the right yields 

j= r w„. = T[w rr \w rr Y:i =1 I4 fc) l 2 



+ ELi E' a =i [/C(^4 fc )*4 fc ) + {if [W„.] - W^ rr • [det 2 sr ]*c«c«* 

= F[w rr ]w rr Y. n k= i l^l'-ELi^l^r-ELi {#W - w^l4 fc) l 2 

+ ELi E:=i [/C(^4 fc )*4 fc ) + {if [W sr ] - • [det z sr } l 2ci k > ■ W„c. 



(B.2) 



(k) 



Using (A. 11), the second term in the last line of R. H. S. of (B.2) is vanished. Substituting 
to K,^) its explicit form obtained from (A. 14), thus, T r W rr is cast into 

7 r W„=7\W rr \W rr YJ l k=l \d 



(fc)|2 



-E 



k=l 



[hN-W rr )^[W rr ]+2{H[W rr ]-E^}] -^ rr |4 fc )| 2 + ELiE: = i/CW4 fc )*4 fc ) (B.3) 

= w rr j[w rr }w rr E n k=1 \4 k ^\ 2 -ELi2 {^[^]-^ (fe) }w- rr |4 fc )| 2 +ELiE: = i/cW4 fc )*4 fc ), 

and taking hermitian conjugate of both sides of (B.3), we have J>W / rr =W / rr J>. This means 
the two hermitian matrices T r and W rr have common eigenvectors to diagonalize them so 
that it leads to (B.l). Therefore, (B.l) and the relation JvW^^W^Jv are equivalent. Due 
to the idempotency relation W 2 r =W rr , J r r W rr =W rr J r r is equivalent to J>lTv r — l^rJvW^^O. 
Further mutiplying (B.3) by W rr from the left and using the explicit form of K$ , we obtain 

W TT T T W rT = W^[W rr }W rr E n k= i I4 fc) | 2 - ELi 2 {H[W rr \ - E^} Wl\c^\ 2 

+ ELi E:=i W rr {(1 2N - W rs )T[W rs ] + H[W rs ] - E^} • W rs • [det z rs ]h c W*cf) (B.4) 

= w rr F[w rr }w rr ELi I4 fc) l 2 - ELi 2 {ffM - iu rr |4 fc )| 2 . 

Subtracting (B.4) from (B.3), it is easy to derive an equivalence relation 

ELi E:=i /CW4 fc )*4 fc ) = F r W rr - W„.T r W„. (B.5) 

Thus, the equivalence of (B.l) with (A. 14) has been proved. The above equivalent relation 
(B.5) also makes a crucial role in the variation of the Res-HB free energy. 
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Let us introduce the following Res-HB free energy which is quite similar to (A. 3) 

but involves the HB interstate density matrix instead of the usual HB density matrix: 

H^Cs (1) = £Li Er, s=1 {H[W rs ] - E^)} ■ [det^4 fe )*4 fe \ 



F[Z] 



HB(2) 
Res 



1 1 

2^ 



T,r, 3 =i Tr i W rs In W rs + (1 2JV - W rs ) ln(l 2JV - W rs )} , 



F[Z]Z {1) + F[Z]^. 



(B.6) 



We are now in a stage to make a variation of the Res-HB free energy. Using the variational 
formulas (A. 12), the thermal Res-HB equation (A. 14), the equivalence relation (B.5) and 
the commutablity [J>, W rr ] = 0, it is made as follows: 



Er=i ^Tr {ELi ESLi K^c^c^uM } 

+ E? =1 ^Tr {^J E2=1 E « =1 K^cf)cf>} 
E"=i ^Tr { [(^W^-W^W^.)^^] + [5u r ul(W rr F r -W rr F r W rr )] } 

~ E"=i Tr [In {H/ rr (l 2iV - W^)- 1 } (1 2JV - W rr )<5lU rr 

+ E' s =l(l2AT - W rs ) ^ {W rs (l 2N - Wrs)- 1 } W rs U r 5u\. + (h.C. 



;b.7) 



(B.8) 



the last of which has no contribution since (1 2 at— W rs )W rs — 0. Then, the variational equation 

5F[Z}™=5F[Z}Z {1) +5F[Z]Z S {2) = leads to In {W rr {l 2N - W^)- 1 } = -/3T r , in which we 

8vl + 5u r ut and 5utu r + vl5u r = 0. Thus we 



Ur 



have used the variational relations 8W r . 
obtain W rr (l 2 N — W^)' 1 = exp{— /3J>} from which we have W rr — exp{— /5J>}(1 2 at — W rr ). 
Finally we can reach the r-th thermal HB density matrix VF rr [S] Res expressed in terms of 

the r-th thermal Res-FB operator J>[Z] Rcs as H / rr [Z] Res = r - — , — - where we 



have used the relation W r 
diagonalized as follows: 

'{Wri} 



^jv+exp^Jv^Res} 
= 0. By using a Bogoliubov transformation g r , W rr [Z] Res is 



W r = r^W rr [Z] Res g r 







1 N - {w ri } 



W r 



-W r 



l+exp{/3e ri } 
1 



,• ••,">), 
(i = l,---,N). 



(B.9) 



l + exp{-/3e ri } 

The diagonalization of the r-th thermal Res-FB operator J>[Z] Rcs by the same Bogoli- 
ubov transformation g r leads us to the eigenvalue e r j. To this eigenvalue by adding a term 
J2k=i{H[W r r[Z} Rcs }—E( k ^ |4 fc )| 2 -l 2A r, the usual HB type of the eigenvalue e ri is realized. 

Using (B.9) we can derive the inverse transformation of (B.9) in the form 

1 



glW rr [J r r \g r = g\ 



W + exp [f3^ r [2] Rcs + U=l {H[W rr [Z] Rcs ]-EW} \Cr k) \ 
{Wri} 



l 2N , 



w, 







■{W ri } 



(B.10) 
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C Calculations of E p A p ,E p B p and J2 P C P at intermediate 
temperature 



First we give a integral formula 



f°° { 1 1 1 7 

/ dy { tanh?/ -— sech 2 ?/ } = — C(3), 

jo [y 5 y z J 7r z 



which can be derived by using the famous mathematical formulas [221 
1 11 1 1 

- tanh y = 8£^ 7 ttt^ ? , -? tanh y - — sech 2 y = 64£°? , 

y y ^ 1 (2m-l) V+Ay 2 ' y 3 U y 2 U ^"^ {(2m-l) 2 7r 2 +4?/ 2 } 



(C.l) 
2 (C2) 



Adopting a new integral variable y— (2m— l)7r/2-tan#, an integral of the second formula in 
(C.2) is easily carried out for hu D 3> 1 as 

64f"^^E- 1 {(2m _ 1) ^ My ^ =32E~ 1 pJ w f^^4c(3), (C.3) 

where we have used £„ =1 (2m-l)- 3 = (7/8) <(3), and C(3) = tt 3 /25.794 36 [22]. 

Next the modified QP energy e is approximated as = \fe 2 + A 2 {2(l±[det^ 12 ]^ /2 )} . 
Let us introduce a new variable y by e = 4^1±[det z^jr^ksTy and quantities xt — A t /?iud 
and y T ±] =^e 2 +A 2 T {4(l±[det z l2 }}]l 2 )y l huj D /k B T where A T = A T N(0)Varcsinh(/iu; D /A T ). 
If £3> At, J2 P \, H v Bp and J2 P C p in (2.14) are recast to the following integrals up to A T : 



2N(0) 



- - ( de 



Jo 



e 2 +A 2 



tanh 



2A; b tJ 



de- 



A 2 T 



V^+A 



: tanh 



\2fc B Ty 



1 (J/^W) 
<1 — ^ 



2^ 



=tanh 



2 +(# } 2t)' 



(C.4) 



'5 A ^ 
de- 



=tanh 



ArE P ^ p 

£_ r 

N(0) ~J-^d £ ^ £ 2 + ^' "~\$k B Tj 2/-R^ e y e2+ A2 



1 /-^d A^ 
cfe- 



=tanh 



2k rT) 



■0, 



(C.5) 



A 2 V (7 /-^o A 2 
2iV(0) Jo 



-r(±) ' 



= tanh 



A/^+A 2 



A1 



T 



<2k B T) Jo eSje 2 +A 



=tanh 



f(±) \ 



2 /-v? 



(y^W) / T dy 



(±)- 



V 



■. tanh 



y 2 + (#^ 



(C.6) 



y 2 \j y^y^xr) 

To get a finite value of X) P A» expanding (C.4) around y^xr, (C.4) is boldly approximated as 

%m^fo T d v\^v-\ (y^T dy tanh ^ sech ^}- (c - 7) 

In a similar way we also get a roughly approximated integral form for (C.6) as 

AtE p Q, /(±)~A 2 r^ J 1 , u 1 
-^(^-(^V)^ dyj^tanhy--: 

Integrations of (C.7) and (C.8) are easily made by using the integration formula (C.l) if we 
take the upper-value y^ to be infinite. 



(C.8) 
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